New ordered phases in a class of generalized XY models 
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It is well known that the 2d XY model exhibits an unusual infinite order phase transition belonging 
to the Kosterlitz-Thouless (KT) universality class. Introduction of a nematic coupling into the XY 
Hamiltonian leads to an additional phase transition in the Ising universality class [lj . In this paper, 
using a combination of extensive Monte Carlo simulations and finite size scaling, we show that 
the higher order harmonics lead to a qualitatively different phase diagram, with additional ordered 
phases originating from the competition between the ferromagnetic and pseudo-nematic couplings. 
The new phase transitions belong to the 2d Potts, Ising, or KT universality classes. 



The low temperature behavior of two dimensional (2d) 
systems with continuous symmetries is controlled by 
topological defects, such as vortices and domain walls. 
Although massless Goldstone excitations, such as spin 
waves, destroy the long-range order of these systems, a 
pseudo-long-range order with algebraically decaying cor- 
relation functions still remains possible. At low tem- 
peratures the topological defects which undermine the 
pseudo-long-range order are all paired up, while above 
the critical temperature these defects unbind, leading to 
exponentially decaying correlation functions and a loss 
of the pseudo-long-range order. A classical example of 
such system is the XY model. At low temperature, the 
topological defects, in the form of integer valued vortices, 
are all joined in vortex-antivortex pairs, resulting in alge- 
braically decaying spin-spin correlation functions. Above 
the Kosterlitz-Thouless (KT) critical temperature 0, Q , 
these pairs unbind and the correlation functions decay 
exponentially. 

Unlike in 3d, for 2d systems the arguments based 
purely on symmetry considerations are not sufficient to 
fix the universality class of possible phase transitions 0- 
3]) and even in 3d a second order phase transition can 
be preempted by a first order one [8|. Violations of 
strong universality are even more common in 2d. Thus, 
it is possible for systems with the same underlying sym- 
metries and the same coarse-grained Landau-Ginzburg- 
Wilson Hamiltonian not to belong to the same universal- 
ity class. It is, therefore, interesting to ask what phase 
transitions are possible for 2d Hamiltonians invariant un- 
der the transformation 6 — > 9 + 2?:. In this paper, we will 
study using extensive Monte Carlo simulations and finite 
size scaling (FSS) analysis, a large class of generalized XY 
models which, while preserving the same — > 9 + 2tt sym- 
metry, have very complex phase diagrams, with phase 
transitions belonging to the Ising and Potts universal- 
ity classes, in addition to the usual KT phase transition. 
In some of these models, transitions can be understood 
in terms of new topological defects, such as fractional 
vortices and domain walls [l|, |9| . Apart from the funda- 
mental considerations regarding the connection between 
symmetry and universality, our purpose is to describe 
new, previously unnoticed, ordered phases which occur 
in 2d systems with continuous symmetry. 

The model considered here has a mixture of ferromag- 



netic and nematic-like interactions, 
H = -^[Acos(^-^) + (l- A)cos(q9 i -q9 j )], (1) 

(ij) 

where the sum is over the nearest neighbors spins on a 
square lattice, < A < 1, and q is a positive integer. 
The first term is the usual ferromagnetic coupling (XY 
model), while the second one favors the adjacent spins to 
have a phase difference of 2fc7r/g, where k < q is an inte- 
ger. Independent of the value of A, the Hamiltonian (TTJ) 
has the symmetry of the pure XY model, recovered when 
A = 1, and is invariant under rotations 9j — > 9j+2ir. For 
A = 0, we have a purely nematic-like Hamiltonian, which 
is also invariant under the transformation 9j — > 9j+2ir/q. 
It is easy to show that in this case there will also be a 
KT phase transition at exactly the same critical tem- 
perature as in the pure XY model. The low temperature 
phase for A = will, therefore, have a pseudo-long-range 
nematic-like order. An interesting question concerns the 
thermodynamics of the model described by Eq. ([1} for 
< A < 1, where both terms compete. 

For q = 2, the Hamiltonian, e q. (ED ), has been stud- 
ied by a number of authors 0, l9rU3j and the presence 
of the second term leads to metastable states in which 
spins have antiparallel orientation. The model has new 
excitations not present in the A = 1 case: half-integer 
vortices connected by strings (domain walls) [l|, across 
which spins are anti-paralelly aligned. At low tempera- 
tures, half-integer vortices are bound in pairs of integer 
vorticity, resulting in a pseudo-long-range ferromagnetic 
order. If A < A mc , as the temperature is raised, the 
string tension between half-integer vortices vanishes and 
the system melts into a nematic phase. On further in- 
crease of temperature, the half-integer vortex-antivortex 
pairs unbind and the system enters a completely disor- 
dered paramagnetic phase. As expected, the transition 
between the nematic and the paramagnetic phases be- 
longs to the KT universality class |l4[. Surprisingly, the 
transition between the two pseudo-long-range ordered 
phases — nematic and ferromagnetic — is found to be 
in the Ising universality class [l[ . This behavior has been 
verified with simulations on the square and triangular 
lattices [ill, E3- On the latter, the geometric frustration 
introduces also a tiny chiral phase above the KT line, but 
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FIG. 1: a) Schematic phase diagram for q — 3. At both extremes A = and 1, the order-disorder transition temperature is 
T c = 0.893. There are two low temperature phases with a pseudo-long range order: generalized-nematic (three preferred spin 
orientations) and ferromagnetic (broken reflection symmetry). Notice the multicritical point around A mc ~ 0.4. For A < A mc , 
there are two transitions: paramagnetic to nematic, in the KT universality class; and nematic to ferromagnetic, in the q = 3 
Potts universality class, b) Schematic phase diagram for q = 8. Besides the paramagnetic, the usual ferromagnetic, and the 
nematic phases, there are two new ferromagnetic phases (Fi and F2) in which spins have half-plane preferred orientations. 



otherwise the phase diagram retains the same topology as 
on the square lattice [12| . A related model, with a similar 
phase diagram, was also studied in Ref. [15j. An inter- 
esting question is whether, for q > 2, the topology of the 
phase diagram remains unchanged. To answer this, we 
have explored using extensive Monte Carlo simulations 
the equilibrium phase diagram of the Hamiltonian ([T]) 
with q ranging from 2 to 10. We find that for q = 2, 3, 4 
the topology of the phase diagram remains the same as 
for q = 2, however, the transition between the pseudo- 
ferromagnetic and pseudo-nematic phases belongs either 
to the 3-state Potts (q = 3) or to the 2d Ising (q = 2, 4) 
universality class, see Fig. [1^. For q > 5, the topology of 
the phase diagram changes completely and new phases 
with a pseudo-long range order come into existence. 

The simulations were performed on a square lat- 
tice of linear size L and periodic boundary condi- 
tions. Both Metropolis single-flip and the Wolff algo- 
rithm [16:] were used. In accordance with the sym- 
metry of the Hamiltonian, the possible order parame- 
ters are mk = L~ 2 |^ exp(ik9i)\ where k = l,.,.,q. 
The corresponding generalized susceptibilities are Xk = 
f}L 2 {{m\) — (mk) 2 ) and the Binder cumulants are Uk = 
1 - (m 4 k )/3(m 2 k ) 2 (H 111. If the transition is not KT, 
the usual FSS can be used to get the critical exponents 
/3,7, and u: to = L~^ v f{tL x l v ) and x = LPi/ v g{tl>/ v ), 
where m is the order parameter and x its susceptibil- 
ity, / and g are the scaling functions, and t — T/T c — 1 
is the reduced temperature. This FSS, however, is not 
valid at the KT transition for which all of the low tem- 
perature phase is critical and the correlation length and 
the susceptibility are infinite Nevertheless, it is 

possible to show that at the KT transition and in the 
low-temperature phase, the critical exponent ratios are 
well defined and the order parameter and the general- 




FIG. 2: Ferromagnetic-nematic transition for q — 3. FSS 
analysis of the susceptibility xi an d the order parameter mi 
(inset) for A = 0.25 and different sizes L. In both cases the 
data collapse is excellent using the 3-state Potts exponents. 
The critical temperature at this point is T c ~ 0.365. 



ized susceptibility scale with the size of the system as 
to oc L~ftl v and x ^ L~*/ v . Exactly at the transition, 
f3/v = 1/8 and 7/1/ = 7/4, which are the same ratios as 
for the 2d Ising model. However, what distinguishes the 
KT transition from the Ising one, is the behavior of the 
order parameter and the susceptibility in the low tem- 
perature phase where they also exhibit FSS, but with 
non-universal critical exponents. Recall that for normal 
second order phase transition, FSS exists only at the crit- 
ical point. This difference can be used to distinguish the 
KT transition from the Ising one. 

FigureQJ, shows the schematic phase diagram for q = 3, 
which is topologically identical to the q = 2 case. At 
high temperatures, the equilibrium state is the disordered 
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FIG. 3: Data collapse for the susceptibility \i an< A the or- 
der parameter 7714 (inset) for g = 8 at A = 0.15 as the sys- 
tem crosses the border between Fi and F2 phases. The col- 
lapse was obtained using the critical exponents of the 2d Ising 
model, 7 = 7/4, v = 1, /3 = 1/8 and T c ~ 0.058. 
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FIG. 4: Angle distribution for the four phases along the fixed 
temperature line T = 0.16 of the q — 8 phase diagram. All 
graphs have the same vertical scale. 



paramagnetic (P). As the temperature is lowered, the 
system enters either the usual ferromagnetic (UF) or the 
generalized-nematic phase (N), depending on the value 
of A. Both of these order-disorder transitions belong to 
the KT universality class. Up to the multicritical point 
located at A mc ~ 0.4, there is a line of critical points 
separating the generalized-nematic from the UF phase. 
The order parameter mi is used to distinguish between 
the generalized-nematic and the ferromagnetic phases: 
mi ~ in the nematic phase and is ~ 1 in the ferromag- 
netic phase. Using FSS, and also the Binder cumulant, 
we find that the critical points along this line belong to 
the 3-state Potts universality class. Fig.[5]shows the data 
collapse of mi (inset) and the corresponding susceptibil- 
ity xii f° r a critical point with A = 0.25. The collapse 
is excellent u sing the critical exponents of the 2d, 3-state 
Potts model [If: /3 = 1/9, 7 = 13/9 and v = 5/6. 

The same topology of the phase diagram persists for 
q = 4, except that the transition from the generalized- 
nematic to the ferromagnetic phase is once again in the 
universality class of the Ising model [2l], [22[. For q > 5, 
however, the topology changes dramatically. Two new 
phases emerge from the T = 0, A = fixed point, see 
Fig. [TJd. The new phases are pseudo-ferromagnetic and 
have a broken reflection symmetry. We shall denote these 
phases as Fi and F2. The low-temperature phase Fi has 
only one preferred spin orientation, while in the phase F2 
there are four preferred spin orientations with different 
weights, see Figs.[T|D and|3J Fig. 2] exhibits histograms of 
spin orientation in the low temperature phases along the 
line T = 0.16 together with a pictorial representation of 
the possible spin orientations for the q = 8 model. All 
the ferromagnetic phases have a quasi-long-range order 
and a broken reflection symmetry. For F2, the distribu- 
tion function has four significant peaks, while for Fi there 
is only a single narrow peak. The UF phase of the XY 
model has a broad continuous distribution of spin ori- 
entations. The nematic phase is characterized by eight 



congruent discrete spin orientations, separated by 7r/4. 
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FIG. 5: Data collapse for the susceptibility xi and the order 
parameter mi (inset) for q = 8 at A = 0.35 as the system 
crosses the border between the generalized-nematic and F2 
phases at T c ~ 0.34. The collapse was obtained using the 2d 
Ising model critical exponents, 7 = 7/4, v = 1 and f3 = 1/8. 

Besides the order-disorder KT transitions, several new 
order-order transition lines appear in the q > 5 phase di- 
agrams, Fig. [TJd. The transition from the Fi to F2 phase 
is well described by the m 4 order parameter, which is 
close to zero in the later phase. The data collapse for 
several lattice sizes is shown in Fig. [3] along with the 
data collapse for the susceptibility Xi- Both collapses 
were obtained using the set of critical exponents of 2d 
Ising model: 7 = 7/4, v = 1 and f3 = 1/8. The tran- 
sition from F2 to nematic is also continuous and is also 
in 2d Ising universality class, as can be seen in Fig. [5] in 
which the data collapse of the order parameter and the 
corresponding susceptibility are shown. For this transi- 
tion, the relevant order parameter is mi. Finally, the 
transition from Fi to UF and the transition from F2 to 
UF are both well described by the ms order parameter. 
The two transitions are found to belong to the KT uni- 
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FIG. 6: Scaling of the order parameter at the KT phase 
boundaries (three bottommost lines) and inside the critical 
phase (topmost line) for q = 8. Notice that along the KT 
transition the exponent /3/v has the Ising value (1/8), while 
inside the critical phase the exponent is non universal. 

versality class with the critical exponent ratios, j/is — 
7/4 and j3/v — 1/8. As an example, Fig. [6] shows that 
F2 is critical with respect to mg order parameter and has 
a non-universal FSS characteristic of a low-temperature 
KT phase, top-most curve of Fig. [5] 

To conclude, we have studied, using extensive numer- 
ical simulations, phase diagrams of a class of general- 
ized XY models described by eq. ([T]). Previous results 
show that for q — 2 besides the usual KT transition, 
there is also a nematic to ferromagnetic transition in the 
Ising universality class. We find that for q = 3 and 4 
the topology of the phase diagram remains unchanged, 
but the ferromagnetic-generalized-nematic transition be- 
longs to the universality class of 3-state Potts and the 



Ising model, respectively. For q = 5, the topology of the 
phase diagram changes dramatically and two new ferro- 
magnetic phases appear. After this, up to q — 10, the 
maximum value explored in this work, the topology of 
the phase diagram remains unchanged. It is very curi- 
ous that although the systems studied here are invariant 
under the continuous global symmetry 9i — > 61 + a (for 
all spins simultaneously), for arbitrary a, the sequence of 
the phase transitions between the pseudo-ordered phases 
follow the one observed for the discrete clock models: 
q=2 Ising; q=3 Potts; q=4 Ising; fo r q =5 a bifurcation 
and a new phase transition appears 23] . If this analogy 
persists, we expect that there should be a critical value 
of q above which the phase transition between Fi and 
F2 becomes KT (24J. This will be explored in a future 
work. The present paper also shows a significant lack 
of universality of 2d systems: while all the Hamiltonians 
studied in this paper have the same underlying symmetry, 
9 — > 9 + 2tt, the transitions between the different phases 
belong to a variety of universality classes. Furthermore, 
since the Hamiltonians discussed in this paper can be 
thought of as the leading orders in a Fourier expansion 
of a general microscopic spin-spin interaction V(9i —0j), 
the work raises a troubling question: How much can we 
really deduce about the thermodynamics of 2d systems 
from the form of their coarse-grained Landau-Ginzburg- 
Wilson (LGW) Hamiltonian? It is clear that the symme- 
try arguments alone are not sufficient to determine the 
phase diagram of these systems, and one needs to have a 
detailed knowledge of the microscopic interactions [7J • 
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